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The band calculation scheme for / electron compounds is developed on the basis of the 
dynamical mean field theory (DMFT) and the LMTO method. The auxiliary impurity prob- 
lem is solved by a method named as NCA/ 2 v', which includes the correct exchange process 
of the f 1 — > f 2 virtual excitation as the vertex correction to the non-crossing approximation 
(NCA) for the f 1 — ► f° fluctuation. This method leads to the correct magnitude of the 
Kondo temperature, Tk, and makes it possible to carry out quantitative DMFT calculation 
including the crystalline field (CF) and the spin-orbit (SO) splitting of the self-energy. The 
magnetic excitation spectra are also calculated to estimate Tk- It is applied to Ce metal and 
CeSb at T = 300 K as the first step. In Ce metal, the hybridization intensity (HI) just below 
the Fermi energy is reduced in the DMFT band. The photo-emission spectra (PES) have a 
conspicuous SO side peak, similar to that of experiments. Tk is estimated to be about 70 K 
in 7-Ce, while to be about 1700 K in a-Ce. In CeSb, the double-peak-like structure of PES 
is reproduced. In addition, Tk which is not so low is obtained because HI is enhanced just 
at the Fermi energy in the DMFT band. 

KEYWORDS: dynamical mean field theory, band theory, Ce metal, Ce-pnictides 



1. Introduction 

The non-empirical band calculation on the basis of the local density approximation (LDA) 
theory is one of the fundamental approaches to obtain the insight into the electronic prop- 
erties of solid state materials. However, it has limitations to treat the dynamical excitation 
spectra (DES) of the strongly correlated electron systems 1 such as / electron compounds. The 
understanding of DES of strongly correlated electrons has been extensively improved recently 
based on the approach of the dynamical mean field theory (DMFT). 2,3 It seems natural to 
develop a non-empirical band calculation scheme for actual materials by the combination of 
LDA and DMFT. 4 - 6 

In the DMFT, the strongly correlated band electron problem is mapped onto the calcu- 
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lation of DES of the auxiliary impurity Anderson model in an effective medium. 3 At present, 
however, we do not have the theoretical method which analytically gives correct DES of the 
impurity Anderson model. Several numerical methods have been developed to calculate DES 
of the Kondo problem. It is known that the quantum Monte Carlo (QMC) method 7 ' 8 and 
the numerical renormalization group (NRG) method 9 give in principle correct DES. They are 
applied to DMFT calculation, 8 ' 10 but they have difficulties in application to the realistic band 
calculation of / electron systems. 6,11 

Approximate but highly flexible approaches have been developed on the basis of the resol- 
vent technique. 12 The NCA method 13 ' 14 has been widely applied to various models though it 
has some weakness in application at very low temperatures: it does not fulfill the Fermi liquid 
relation. But it can be easily applied to realistic situations such as the competition between 
the Kondo effect and the crystalline field (CF) splitting. 1 NCA has been also applied to the 
DMFT calculation. 15 The NCA equation is justified by the 1/N f expansion when the fluctu- 
ation of the valence is restricted to the f 1 and f° configurations, where Nf is the degeneracy 
factor of the / orbital. 14 However, the exchange coupling through the virtual excitation to 
the f 2 configuration is not negligible in quantitative calculation. A scheme which includes the 
f 2 configuration in the frame work of the NCA type diagram has been used. However, it does 
not lead to the Kondo temperature (Tk) given by using the exchange constant obtained by 
the Schrieffer- Wolff (S-W) transformation. 16 Therefore, the exchange process through the f 2 
configuration is not properly accounted in the scheme. Actually, the DMFT calculation by 
using this method gives a too small Kondo resonance peak compared with the result using 
QMC method. 6 In the DMFT band calculation, it is desirable to use a method which gives 
the correct order of magnitude of 7k . 

A method to include the f 2 configuration within the 1/Nf order has been developed. 17-19 
The NCA equation is modified by a vertex correction term which leads to exchange cou- 
pling through the f 2 configuration. The self-consistency equation of this method needs huge 
computational time, if one tries to solve it strictly. But we can usually apply simplifying ap- 
proximations for the equation. 17 It does not need so much computational time, only few times 
larger than that of the NCA calculation, and gives Ik consistent with the result of the S-W 
transformation. We call the method as NCA/ 2 v', henceforth, and is explained in Appendix. 

Several efficient band calculation techniques have been developed. 20 Among them the 
LMTO method is a very convenient one because it has formal similarity to the usual LCAO 
picture. 21 It is widely used though accuracy is not so very high when it is compared with 
other sophisticated methods. In this paper we develop a DMFT band calculation scheme by 
the combination of NCA/ 2 v' and LMTO methods, and apply it to Ce metal and CeSb as the 
first step. 

The effective hybridization intensity (HI) is defined nominally as the product of the square 
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of the hybridization matrix and the density of states of conduction bands in the auxiliary 
impurity problem. It plays the key roles in the DMFT calculation. We mainly study how HI 
is modified in the DMFT calculation from that of LDA theory, and how the effects of the 
spin-orbit (SO) splitting and the CF splitting appear. 

The relative occupation number of excited SO level is strongly reduced from that of LDA 
in the correlated bands even when the total occupation number of the 4/ electron is fixed 
near the value of LDA result. 

Ce metal is known to have iso-structural fee phases. 22 The low temperature state is called 
as the a phase and it has smaller lattice constant. The higher temperature one is called as 
the 7 phase. The DMFT studies on this transition have been extensively carried out. 6 ' 15 ' 23 ' 24 
In these studies the SO and the CF splittings were neglected. In the present study we are not 
concerned with the a — 7 transition, but we calculate the photo-emission spectra (PES) 25-28 
for the fee state with the lattice constant of these two cases. The HI estimated by the LDA 
band has a peak just below Ep. This peak is reduced in the DMFT band, thus the SO side 
peak in PES becomes relatively conspicuous. 

The CF level splitting is small in the auxiliary impurity model of 7-Ce, and Tk is estimated 
to be about 70 K by the calculation of the magnetic excitation. In a-Ce, Tk of the auxiliary 
impurity model becomes very high about 1700 K. 

CeSb would be a first material for which HI calculated by the realistic band is used in 
the analysis of PES. 26 ' 29 ' 30 The hybridization matrix is not small in the valence band region 
of Sb p states. In the energy region near Ep the density of states of bands is small reflecting 
the semi-metallic nature of the band structures. Various anomalous properties such as the 
very large magnetic anisotropy and the origin of the ferromagnetic layers in the ordered state 
are explained on the basis of the p — f mixing model proposed by Kasuya et. a/. 31-34 The 
top of the valence p band has the (j = 3/2)Fs symmetry. A part of T$ bands which has the 
same symmetry with that of the occupied 4f band is pushed up above Ep due to the p — f 
hybridization in the ordered phase, and this plays important roles in the p — f mixing model. 
Such re-constructed electronic band structures in the ordered states are confirmed by the 
dHvA and the optical experiments. 35-37 In the ordered states we can use the Hatree-Fock- 
like static mean field approaches, and detailed comparison with experiments has been carried 
out. 34,38 However, in the paramagnetic phase we need DMFT band calculation. 

CeSb has a characteristic double peak structure in PES; the shallow-energy one is at 
about 0.8 eV below Ep and the deep-energy one is at about 3 eV below £? F . 25 > 30 > 39 > 40 This 
has been explained by the large peak of HI in the valence band region. 26,29 ' 30 Tk has been 
calculated by using the same HI on the basis of the impurity model. 41 ' 42 It is very low, about 
10~ 5 K because of the small density of states near Ep. However higher Tk, about 10 K, has 
been expected from the transport properties. 32 In the present DMFT calculation, the double- 
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peak- like structure is obtained. At the same time we obtained a rather higher Tk, since HI 
of DMFT has a sharp peak just at Ep caused by the hybridization of the bands with the 
correlated 4/ bands. Similar evolution of the peak at Ep in HI has been observed recently 
in DMFT calculation for CeP and CeAs in ref. 43. But the SO splitting and the correct 
exchange process through the f 2 state have not been included there. The (j = 5/2)rs state 
has lower energy than that of (j = 5/2)Tj in the present calculation starting directly from the 
LDA band. It is demonstrated that more similar result with experiments such as the definite 
double peak structure of PES and lower, but not so low, Tk can be obtained when the band 
overlapping between the occupied Sb p bands and the un-occupied Ce d bands is reduced from 
that of the LDA band. 

In both materials, we find that HI which is directly calculated from the 4/ partial DOS of 
LDA bands seems to be somewhat large when we compare the calculated PES with those of 
experimental ones. 25 ' 27 ' 41 But in the DMFT calculation the effect of HI is relatively weakened. 
The calculated HI is somewhat lager, but not so drastically large than the magnitude expected 
from experimental data. 

In §2, we give our formulation on the basis of the LMTO method. The DES which are 
calculated by the NCA/ 2 v' method for the single impurity model are compared with the 
results of NRG in §2. Applications to Ce metal with the lattice constant of 7-Ce, and with 
that of a-Ce are given in §3 and §4. In §4 the calculation for CeSb is shown. Summary is given 
in §6. In Appendix we explained the NCA/ 2 v' method. 

2. Formulation 

2. 1 LMTO-DMFT Matrix Equation 

We consider the spectra of the following Hamiltonian, 6 

n = n LDA + u - £(£ <£. Tl c ri r7 - nf/) 2 . (i) 

i r, 7 

Here r ^ is the annihilation operator for the atomic localized state, 0|p 7 ( r )> at site i 
with the 7 orbital of the T irreducible representation. The quantity = ^2r-y n ]^ip~ i * s 

the occupation number on the atomic 4_f -electron per a Ce ion in the LDA calculation. We 
assume that the local Coulomb interaction is expressed by the (f>^ p ^ orbit. 
The excitation spectra are expressed by introducing the self energy terms, 3 

Hdmft = Hlda + Yl ^(e + iS) + ep - ^ DA )|(A| r7 >< <t>\ r J, (2) 
i ,(r l7 ) 

where e r is the single electron energy level of 4/ state, and £p DA is the energy level in the 
LDA calculation. The self-energy Sr(e + id) is calculated by solving the auxiliary impurity 
problem in the effective medium with the use of the NCA/ 2 v' method, 17 the outline of it is 
described in Appendix. 
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Usually, the DMFT calculation is carried out after transforming the TihDA part into the 
LCAO scheme or to the localized Wannier representation, and then using those basis func- 
tions. 6 In this paper we include the self-energy term directly in the LMTO matrix. This 
simplifies the calculation because the processes of defining the localized basis function are 
skipped, but instead, Greenian equation in the non-orthogonalized bases must be introduced. 
In later calculation we will approximate the localized 4/ state a by the band center orbit 
<p(—) which has logarithmic derivative — £ — 1 on the muffin-tin surface, because it it well 
localized for the 4/ state. 

In the LMTO method, 21 the Hamiltonian Wlda is represented as the matrix using the 
LMTO bases, 

^ fc (r) = XXSx?L(r). (3) 

qL 

ik 

Here, a L is the expansion coefficient of the j-th eigen state with the wave number k on the 
LMTO base of the angular momentum (I, m) and spin {a) at site q in the unit cell. We denote 
as L = (£,m,a), henceforth. The LMTO Bloch state is given as 

Y k (r) _ * q U-,r q ) v 3Vl'(+;*V) k U) 

x " L[r j " ^V2<M-) 99 V (2* + i)y/stf2*M+) q ' L '^ 

in the sphere of q', where r q < = r — R q < is the electron coordinate from the center of MT at 
R q i. The quantity S q is the radius of the MT sphere, and S^ L , qL is the structure factor. 21 
The LMTO base function with logarithmic derivative D is given 

% L (D; r ) = i e Y e m (r ,) X (a)(^(r 9 ) + w(D)^(r,)), (5) 

with the omega factor oo(D) of Andersen. 20 The quantity 3> 9 l( — ; r q ) (<& q i{+; r q )) in eq. (4) 
means (^ + l);r q ) ($ q L(£;r q )), and $^(±) is (f> u e(r q ) + uj(±)<j> u e(r q ) on the MT sphere. 

The wave function (j) v i{r q ) and its dot state, 4> u e(r q ) are calculated at the energy e v . 

The density of states (DOS) of the excitation is obtained by solving the Greenian equation 
of the matrix form, 

[z6 - Hlda - t{z)]G{z) = I. (6) 

The notation hat of the matrices means that they are defined on the bases of x^l- ^ s noted 
previously, the LMTO bases are not orthogonalized with each others. 21 Therefore, we have 
the overlapping integral O in eq. (6), whose matrix element is given as 0^ L , qL = (x^l' I ^qi)- 
The DOS on the /-state is given by 

/0 | band )( £ ) = -Itr[O r G(e + i<5)]. (7) 



7T 



The projection operator is defined as 



Or = E^ ( r 7 )><^ ( r 7 )l- ( 8 ) 
i 1 
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The DMFT self-consistency calculation is carried out by the following way. (I) At first we 
solve the auxiliary impurity model with a trial HI and impurity level e r , and calculate the 
DOS /9p np '(e), then calculate the local impurity Green's function by the Cauchy integral, 

1 „ (imp.)/ \ 

G^\z) = - / dx^ ^, (9) 

gr J z-x 

where gr is the degeneracy factor of the V symmetry state. (II) Then the electron self-energy 
due to the Coulomb interaction, Er(z) is calculated by the relation, 

E r (*) = z - 4 imp ° - G^iz)' 1 - xP(z). (10) 

Here Ep 1 ^ (z) is the hybridization self-energy of the impurity problem, which is calculated by 
the Cauchy integral from the effective trial HI. 44 (Ill) By using the self-energy Sr(^) in (6), 
the band problem is solved, and using the equation (7), the DOS of the band p r band ^(£) is 
calculated. (IV) If p r bEmd) (e) and p r imp °(e) agree with each other within the tolerance, the 
calculation is stopped. (V) If both do not agree, local band Green's function, G r band \z) is 
calculated by the Cauchy integral, and we estimate new effective hybridization self-energy 
from the relation, sj! 1 ' 116 ™^.?) = z — e^ S '^ — G ( ^' and \z)~ 1 — Er(z), where e^ S '^ is the first 
moment of p r ba ' nd \e). (VI) Then new effective HI is calculated from the imaginary part of 

s (h:ncw)^ + We retum tQ gtep ^ 

We note that the spectrum intensity in the resolvent method is not usually normalized to 
be unity, because configurations of the 4/ state are restricted to some important ones. 14 This 
is overcome by introducing the normalization factor and the level shift Aer as 

E r (z) = Ae r + Sr(*) - (Zf 1 - l){z - Bp), (11) 

where is the total integrated intensity of DES, and Aer is determined from the first 

moment of p^ mp '\e). These ensure the relation Er(-z) ~ 0(1/ z) for large \z\, which is re- 
quired from the usual Kramers-Kronig relation. This means that the eq. (9) is expressed as 
G^\z) = [z- 4 imp -) - Ae r - E r (*) + ^ - l)(z - E F ) - E^(z)]- 1 . 
2.2 NCA f 2 v' Method 

Before going to the DMFT calculation, we compare briefly the solutions to the impurity 
problem obtained by four typical methods of resolvent technique. 17 By the solid line in Fig. 1, 
we show the single particle spectra calculated by the NCA/ 2 v' method for simplified hybridiza- 
tion. We have assumed that: the conduction band has constant density of states between the 
energy region (-D ~ D) with D = 5; HI is v 2 p c = 0.08 with the hybridization matrix v and 
the density of states of the conduction band p c ; the doubly degenerate /-states are located 
at e/i = —2.2 and e/2 = —1-9; the Coulomb interaction constant between the /-electrons is 
U = 6. The temperature is set at T = 1.723 x 10~ 3 . There appears the Kondo resonance peak 
just at the Fermi energy, and satellites on both sides of it. The satellites correspond to the 
shake up excitation e/i — > e/2 in the creation of hole or electron near the Fermi energy. The 
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result obtained by usual NCA for the (J -/ 1 ) fluctuation model is given by the dashed line. 

The Kondo resonance peak is small because the exchange coupling through the f 2 configu- 

i 

ration is not considered. Grossly Tk is estimated as De N f Jpc with the exchange interaction 
constant, J. The constant NfJp c is given as NfV 2 p c /\ef\ in the usual NCA, where |e/| is the 
excitation energy of f 1 — > f°. When the virtual excitation to the f 2 states is included in the 
S-W transformation, 16 this factor is replaced by the expression N fv 2 p c /\ej\ + NfV 2 p c /\ej + U\. 
The quantity \ef + U\ is the excitation energy of f 1 — > f 2 . Though we usually expect the re- 
lation |e/| < \ef + U\, the second factor gives a comparable magnitude with the first term. 
Because the inverse of J appears in the exponential function, the small increase of J greatly 
increases Tk. 




Fig. 1. Comparison of spectra calculated by several methods on the basis of the resolvent technique. 
The dashed line is the result of usual NCA. The two-dots-dashed line is given by NCA type 
diagram including the / 2 states. The solid line is given by NCA/ 2 v' and the dot-dashed line is 
given by NCA/ 2 v. Two doublets are located at e/i = —2.2 and e/2 = —1.9. The Coulomb constant 
is U = 6 and the hybridization intensity is v 2 p c — 0.08, the density of states of band is constant 
between (-D ~ D) with D = 5. T = 1.723 x 10~ 3 . The DOS of this figure is normalized as the 
states/spin. The zero of the lines are shifted by 0.2 successively, but the dot-dashed line by 0.4 
from that of the solid line. 



If we tentatively use Nj = 2 by assuming the energy difference e/2 - e/i to be large, Tk 
is estimated as Tk,s-w = 7.0 x 10~ 4 and Tk,nca = 5.3 x 10~ 6 . If we use Nf = 4 by assuming 
that the energy difference is small, Tk is estimated as Tk,s-W = 0.073 and Tk,nca = 0.0085. 
These estimated values are largely different with each other. In quantitative calculations, it 
is important to calculate the DES by including the f 2 configuration, and also the SO and CF 
energy splitting. 

The spectra obtained by the sum of the NCA type diagrams for the (Z -/ 1 -/ 2 ) fluctuation 
model is given by the two-dots-dashed line. The Kondo resonance peak in this method is 
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smaller than that of the usual (Z -/ 1 ) type NCA calculation. We think this result is caused 
by the relative decrease in J 1 state energy: the virtual excitation to the f 2 state has order 
Nj — 1 effect, and thus decreases energy in f 1 state largely; this decreases the exchange 
coupling of the virtual excitation to the /° state. On the other hand the exchange coupling 
through the virtual excitation to the f 2 state is not included in this approximation. Therefore 
Tk decreases. 

We show by the dot-dashed line the spectra calculated by a method named as NCA/ 2 v. 
Details of which are explained in the Appendix. In this method, the off-diagonal contribution 
of the vertex is added to the spectra of NCA/ 2 v'. It gives the interference effect of /° <-> f 1 
and f 1 <-> f 2 processes in the final state of the single particle excited states. This effect 
transfers some intensity from the high excitation energy region to the low energy region, but 
total intensity is not changed. 17 The correction is not so drastic when compared with the 
change from NCA to NCA/ 2 v'. The calculation of this off-diagonal correction needs the very 
large computing time of order N £ , while the time for diagonal part needs only N 2 , where N £ 
is the number of mesh points of energy. N £ ~ 10 5 is needed because the width of the energy 
mesh should be smaller than the Kondo temperature. We use the NCA/ 2 v' method in the 
DMFT band calculation in this paper. 




Fig. 2. Comparison of spectra calculated by NC A/ 2 v' (dashed line) and NRG (solid line) methods. 
The temperatures are T = 0.7 x 1CT 4 in NCA/V and T = in NRG. The parameters are the 
same as those in Fig. 1. The discretization parameter of NRG is A = 2. The peak height at the 
Fermi energy of NCA/ 2 v' method is about 1.1. 



In Fig. 2, we compare the spectra obtained by NRG method 9 (solid line) at T = and 
the NCA/V method (dashed line) at T = 0.7 x 1CT 4 . This temperature is about 1/10 of 
the Kondo temperature, Tk,sw(Nj = 2) = 7 x 10~ 4 , and will be low enough to compare the 
result of NRG at T = 0. The intensities of the peaks at the Fermi energy are comparable, 
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slightly small in the NCA/ 2 v' method. The integrated intensity of the shake up side peaks 
is considerably small in the NCA/ 2 v'. This points will be improved if we use the NCA/ 2 v 
method, but we restrict ourselves to the NCA/ 2 v' in this paper. We note, even the intensity of 
the shake up peaks is small, the present method gives the energy scale of Kondo effect not so 
different from the calculation of NRG. Usual NCA and NCA/ 2 give some order of magnitude 
lower Kondo temperature. 

3. 7-Ce 

3.1 The hybridization intensity of LDA band 

As noted previously we approximate <^>| r ^ by <^ r ^(— ), and e v = e r LDA ^ is chosen to fulfill 
the condition u>(— ) = 0. The energy e r = is changed as a parameter to lead the target of 
the total 4f occupation number. In Fig. 3, we show DOS calculated by the LDA band for the 
lattice constant of 7-Ce (a = 9.75477 AU). The 5p of Ce is treated as the frozen core state and 
the scalar relativistic approximation is used. The mesh of k-points in the band calculation is 
generated by taking 9 points on the T — X axis (including T and X points). The SO interaction 
is added as (es ■ I in the matrix. The solid line in the top panel gives the total DOS. The 
shaded region in the panel gives the partial DOS on the (j = 5/2)IV, the middle panel gives 
the partial DOS on the (j = 5/2)Fg, and the bottom panel gives partial DOS on the j = 7/2 
state. The Fermi energy Ep = 0.4585 Ry is indicated by the vertical solid line. For the j = 7/2 
group, DOS is summed over CF states because their splitting is not so important since they 
have higher excitation energy in the DMFT band. In the inset, the partial DOS's near Ep 
are shown for the (j = 5/2)r7(dashed line), (j = 5/2)rg(dot-dashed line) and j = 7/2(dotted 
line) states. 

The total occupation number of the 4f electron is obtained as 1.0859, and partial occu- 
pation numbers are shown in Table I. Hereafter we denote ^ p 7 (— ) as 4>^ r ^. By using DOS 
Pp LDA ^(e) on the orbit 4>{ rj , we define Green's function G r LDA \z) by the Cauchy integral 
same to eq. (9). Following ref. 44, the hybridization self-energy for the non- interacting case is 
given as Ep h ' LDA ' ) (z) = z — e |f ff ' LDA ) _ G r hDA \z)~ 1 . The HI is obtained from the imaginary 
part of 4 h ' LDA) (£ + W)- 

The HI for the LDA band is shown in Fig. 4. We note that HI of (j = 5/2)T 7 (upper 
panel) is relatively large near Ep where the 5d band has the large density of states, while it 
is small in the energy region of the deep valence band where the 6s band is dominant. In the 
inset of Fig. 3, we have shown DOS for /-components near Ep. The DOS of (j = 5/2)T7 is 
relatively larger just below the Fermi energy than that of the (j = 5/2)r§. The degeneracy 
factor of the state is half of the (j = 5/2)rg. Therefore, HI of the (j = 5/2)r-j is about twice 
of the (j = 5/2)r§ component even when DOS has comparable value. This fact can be seen 
in Fig. 4. The HI of the j = 7/2 has the comparable magnitude with that of the (j = 5/2)T7 
near Ep, and with that of (j = 5/2)rg in the deep region. The averaged HI over these three 
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Fig. 3. The density of states (DOS) obtained by LDA band calculation for 7-Ce. The solid line in 
the top panel is the total DOS, the shaded region is DOS on the (j = 5/2)r 7 , the middle panel is 
DOS on the (j = 5/2)r 8 , and the bottom panel is DOS on the j = 7/2 states. The Fermi energy 
Ep = 0.4585 Ry is indicated by the vertical solid line. In the inset, DOS's near Ep are shown 
for (j = 5/2)r 7 (dashcd line), (j = 5/2)r 8 (dot-dashed line) and j = 7/2(dottcd line) states. The 
standard LMTO method with ASA and CC described in ref. 21 is used. The exchange correlation 
function of Gunnarsson-Lundqvist 45 is employed and the lattice constant is a = 9.75477 AU. The 
frozen core approximation is used for 5p states. 



groups is similar to that used in ref. 28. 

When HI is strong, the /°-like peak in DOS of the correlated band appears below the 
band bottom (0.2020 Ry). In such case the peak becomes very sharp, thus it is not easy to 
proceed the iteration steps. To avoid such situation, we added small HI in the energy region 
below the band bottom in the self-consistent DMFT calculation. 46 The added part is included 
in Fig. 4. 

The 4/ spectra calculated by using the NCA/V method at T = 300 K for this hy- 
bridization is shown in Fig. 5. The levels are chosen to be equal to those obtained in the 
self-consistent calculation. The solid line gives the total 4f PES. The single particle DOS's 
for each component is also shown. Because the multiplet splitting of 4/ 2 states is neglected, 
we have only single large peak at about Sf + U ~ 0.8 Ry for the f 1 — > f 2 excitation. The 
Coulomb interaction constant is assumed to be U = 6.5 eV (= 0.48 Ry). 

The occupation numbers are shown in the row of nj! mp in the Table I. The occupation 
number of j = 7/2 components is reduced to 0.0949 from the value, 0.3997 of the LDA 
calculation. Even when the total occupation number of the 4/ electron is not so different from 
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Fig. 4. Hybridization intensity (HI) calculated by the LDA band for 7-Ce. The upper panel is HI for 
(j — 5/2)Ti, the middle panel is one for (j — 5/2)r§, and the lower panel is one for j = 7/2 state. 
In the inset, Hi's near E-p are shown for (j — 5/2)r7(dashed line), (j — 5/2)Pg(dot-dashed line) 
and j = 7/2(dotted line) states. Small HI is added below the band bottom (at 0.2029 Ry), see the 
text. 



the value of LDA, the partial occupation will be largely changed in the strongly correlated 
bands. 

The relative occupation of the (J = 5/2)Ti states to the (j = 5/2)r§, 0.672, is large com- 
pared with 0.341 which is expected from the thermal population with their energy difference, 
115K. Larger HI near Ep enhances the population of (j = 5/2)T7 when the HI of the LDA 
band is used. 
3.2 DMFT calculation 

In Fig. 6, we show DOS of the 4/ components obtained by the DMFT band calculation 
at T = 300 K. 

The effective HI is shown in Fig. 7. Compared with the result in Fig. 4, the peaks of HI 
at about 0.015 Ry below Ep are reduced. At the same time the peaks of (j = 5/2)r§ and 
j = 7/2 components are pushed up just above E-p. 

The HI in the DMFT has the small peaks (or shoulder for the j = 7/2 component) just 
at Ep. The averaged HI near Ep is reduced. Reflecting this the peak of the spectra at Ep is 
relatively weak, and the SO side peak in PES becomes conspicuous in the DMFT band. Such 
reduction of HI just at Ep has been also observed in previous calculation. 15 

We started calculation by choosing the energy levels of (j = 5/2)T7 and (j = 5/2)rg to 
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Fig. 5. 4/ spectra of the impurity model by using HI obtained by the LDA band for 7-Ce at T = 300 
K. The solid line is the total 4/ PES spectra. The dashed line is DOS of (j = 5/2)T 7 , the dot- 
dashed line is DOS of (j = 5/2)P 8 and the dotted-line is DOS of j = 7/2 components. Inset shows 
the spectra near the Fermi energy. 



Table I. Various quantities obtained in DMFT calculation for 7-Ce. rip DA is the occupation number 
in the LDA band, rip lmp '•' is the occupation number obtained by the impurity problem using HI 
of the LDA band. n r band ^ is the occupation number in the DMFT band. 4 imp ' } is the occupation 
number in the auxiliary impurity problem of effective HI of DMFT calculation, er is the energy 
level, Zr is the renormalization factor of the / band, £r is the effective energy of the renormalized 
band, and Tr is the imaginary part of the self-energy. The energy levels are measured from the 
Fermi energy E F = 0.4585 Ry, and given in Ry. e\, = -0.12170 Ry. (4/ = 0.007032 Ry. 
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be equal, -0.1500 Ry. The energy level of j = 7/2 is initially set as -0.1300 Ry. At the first 
stage, (A) HI is fixed as HI of LDA, and only energy levels are changed. We search for the 
solution which gives the partial occupation numbers of band with the self energy grossly agree 
with those of the impurity model. The target of the total occupation number of 4/-electron 
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Fig. 6. 4/ spectra of the auxiliary impurity problem for 7-Ce. For the definitions of lines, see the 
caption of Fig. 5. 




Fig. 7. Effective HI in DMFT for 7-Ce. For the definitions of lines, see the caption of Fig. 4. 



is set as 0.982. In the next stage, (B) the self-consistent spectra and the occupation numbers 
are determined by changing HI and levels by the iteration method. We obtain energy levels 
in the auxiliary impurity model, E((j = 5/2)r 7 ) = -0.11943 Ry, E((j = 5/2)r 8 ) = -0.12016 
Ry and E(j = 7/5) = —0.09595 Ry. Even when we skip the stage A, we can obtain similar 
self-consistent results, but not exactly identical ones in the strict sense because the tolerance 
is not small. 

The occupation numbers in the auxiliary impurity problem and the DMFT band are 
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shown in the Table I. The tolerance for the relative difference of 4/ DOS is set to be 0.1 for 
the j = 5/2 group and 1.5 x 0.1 for the j = 7/2 component. 

The energy level of j = 7/2 states is about 0.32 eV (= -0.09595 - (-0.11943) Ry) higher 
than the energy of j = 5/2. The SO splitting is slightly enhanced than the value, 0.29 eV 
estimated by the SO constant in the LDA band calculation. 

We obtained the cubic crystalline field splitting in the auxiliary impurity model, 115 K= 
( = —0.11943 — (—0.12016)) Ry, which is caused purely by the hybridization effect because we 
have not included the electrostatic CF field in the calculation. However, this difference should 
not be considered literally. The difference in HI will also cause the effective CF energy. The 
relative occupation number of (j = 5/2)T7 to (j = 5/2)rs, 0.3292/0.6165 = 0.53 is almost 
equal to the ratio of their degeneracy factor (0.5), but somewhat larger than it. This means 
that their effective energy levels are almost equal, but (j = 5/2)T7 has slightly lower energy. 



1 . 50 




logioZs 

Fig. 8. Magnetic excitation spectra calculated by using the effective HI for various compounds. The 
horizontal axis is the logarithm of energy in Ry, and the spectral intensities are normalized by 
the values at the lowest energy (E = 10~ 7 Ry). The bold (thin) solid line is calculated for 7-Ce 
(a-Ce). The bold dot-dashed line is calculated for CeSb. The bold two-dots-dashed is calculated 
for CeSb with reduced Sb p - Ce d band overlapping, see the text. For these lines, the full matrix 
elements of magnetic moment in the j = 5/2 manifold are used. The bold dashed line is calculated 
for the fictitious model of 7-Ce in which the matrix elements of magnetic moment within the 
(j = 5/2)T 7 manifold are used. The thin dot-dashed line is calculated for the model of CeSb with 
matrix elements within the (j = 5/2)r 8 manifold. The thin two-dots-dashed line is calculated 
for the model of CeSb with the reduced band overlapping and the matrix elements within the 
(j = 5/2)r 7 manifold. 

In Fig. 8 we show the magnetic excitation spectra of the auxiliary impurity model, 
^9x(a;) 1 _ c - / 3^ , as a function of logarithm of energy (in Ry). The intensity is normalized 
by the value in the low energy limit. The bold solid line and the bold dashed line are calcu- 
lated for the effective HI for 7-Ce at T = 300K (i.e., Fig. 7). The bold solid line is obtained 



14/30 



J. Phys. Soc. Jpn. Full Paper 

by using the matrix element of the magnetic moment within the j = 5/2 manifold, while the 
bold dashed line is calculated by the fictitious model in which we use the matrix within the 
Tj manifold. The bold solid line decreases to half of the low energy value at —3.37, which 
corresponds to the excitation energy about 67 K. This energy will be ascribed to the measure 
of Tk- The bold dashed line decreases to half of the low energy value at slightly higher energy 
side, -3.32 (76K), but the energy difference is very small. When we consider the fictitious 
model with the matrix elements only within the T$ states, the magnetic excitation spectra 
shift to low energy side by about 10 K. The CF splitting energy is small compared with Tk- 

We have also calculated the magnetic excitation for the case of Fig. 5. It decreases to half 
of the low energy limit at —2.6 which corresponds to 400 K. Tk is reduced in the DMFT band. 

The mass renormalization factors and the effective levels are shown in the Table I. Here 
we have approximated as e r imp '' > + Aer + Sr(e) — (Z^ 1 — l)(e — Ep) ~ er — (Zp~ 1 — — Ep). 
The effective levels are measured from the Fermi energy. We found DOS by using these 
renormalized band parameters shows, only in the restricted energy region very near Ep, similar 
value to DOS obtained by the direct DMFT calculation. The renormalization factor about 40 
for j = 5/2 states is very large. But we think these values should not be considered literally 
because the temperature T = 300 K is very high compared with the Kondo temperature. The 
imaginary part of the self-energy at the Fermi energy is estimated to be about 2 eV. The 
actual effect will be renormalized as 1/40 ~ 0.05 eV, but is still not small. 

The change of the band structure in DMFT will also be important in the quantitative 
calculation of the specific heat. The 4/ component at the Fermi energy in the DMFT band 
(about 6 states/Ry, see Fig. 6) is small compared with that (about 35 states/Ry, see Fig. 3) 
in LDA. If one define the mass enhancement factor as the ratio to the DOS of the LDA 
calculation, the effective factor will be reduced by the quantity, 

At this point we note the difficulty in DMFT calculation. The self-energy of the hybridiza- 
tion is relatively very small, about 1/100 of the self-energy due to the Coulomb interaction 
except the energy region near Ep. The latter is roughly estimated as Z~ l NfV 2 p c in the resol- 
vent method, and we have Z~ 1 Nf ~ 500, where Z^ 1 is the mass enhancement factor. This 
causes trouble in the iteration procedure of DMFT calculation. Usually the new trial HI is 
estimated from the difference of the quantities Gp 1 (e) = e — ep — Sp -* (e) — T,^ (e) obtained by 

the band calculation and the impurity problem for the trial HI. (In this place, the self-energy 

(c) 

due to the Coulomb interaction is denoted as £ r '(e) instead of Sr(e) to avoid confusion.) 
Usually we have relation |£ r C ^(e)| >> |Sp h ^(e)|, then the very small uncertainty of E r C ^(e) 
greatly affects new trial of S r h ^(e), while the self-energy Sp ^) itself is very sensitive to the 
detail of HI. These facts lead to the unstable iteration loop. The change of the new trial HI 
from that of the preceding step must be restricted to be very small. Therefore, the convergence 
to small tolerance is very slow. The fact |x£ C) (e)| » |£ r h) (e)| means, in other words, that 
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any trial HI usually gives the almost self-consistent solution except very near Ep. 

We do not allow the large change of HI except for the energy region near Ep. We also 
stopped the DMFT iteration steps less than 60 though the maximum relative difference of the 
spectral density is not so small, about 0.1 after the smoothening process of spectra. Usually 
difference appears in the steeply increasing (or decreasing) region of the sharp peaks of the 
spectral density. The difference can not be reduced easily by further iteration steps. However, 
the difference of the peak height itself is usually small. In addition, the relative difference 
of the densities just at Ep usually becomes less than 0.01. The tolerance for the j = 7/2 
component is set about 3 times larger because maximum difference of it appears near the SO 
side peak in the unoccupied energy region not near Ep . 

To give the atomic 4/ electron number to be near unity (0.983), we obtained the 4/ level 
to be about —0.120 Ry below the Fermi energy. The calculated 4/ peak locates at about 
2.2 eV below Ep, which is slightly deep compared with the energy 2 eV observed in PES 
experiments. 25 ' 27 ' 28 If we need the 4/ level to appear near the experimental position, the 4/ 
occupation number must be chosen to be 0.97. 

In preceding DMFT calculations, a hump at Ep is obtained by QMC calculation, 6 but 
the sharp structure has been smeared out because of the broadening and higher temperature 
calculations. The Kondo peak has not been observed in NCA calculation. In experimental 
studies, SO side peak similar to the result in Fig. 7 is observed in the high resolution ex- 
periments. 25 ' 27 ' 28 In the present calculation we have obtained that the CF splitting energy is 
smaller than Tk. But this result should be checked carefully because we have neglected the 
electrostatic energy part of the CF splitting. 

Detailed comparison with experiments will be carried out in the near future. 

4. a-Ce 

Next we study the Ce metal with the lattice constant of a phase (a = 9.16517 AU). 
The total occupation number of the 4/ electron in the LDA is 1.053 which is almost equal 
to that of 7-Ce, but is 0.033 smaller. We choose the target of the 4/ electron number as 
0.95 = 0.983 — 0.033 in a-Ce. The obtained occupation numbers and the energy levels are 
summarized in Table II. The 4/ DOS is shown in Fig. 9. The small SO side peak is identified 
in PES. This result is similar to the result in experiments, but the position of the f° excitation 
peak is very deep compared with experimental results. 

The magnetic excitation spectra are shown by the thin solid line in Fig. 8, which show 
peak at —2.7 (310 K) and decrease at —1.98 (1650 K) to half of the low energy limit. When we 
use the fictitious model with matrix elements within the Tg manifold, the spectrum intensity 
of the peak region is reduced somewhat. 

In Fig. 10, we show the effective HI. Peaks of HI near Ep in LDA are reduced to the 
narrow peaks just on Ep in the DMFT band. Thus Tk becomes lower than that obtained by 
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Fig. 9. 4/ spectra of the auxiliary impurity model for a-Ce. For the definitions of lines, see the 
caption of Fig. 5. 



Table II. Various quantities obtained in DMFT calculation for a-Ce. Ep — 0.5796 Ry, and energy 
levels are given in Ry. efy = —0.16200 Ry. £4/ = 0.007058 Ry. For the definition of quantities, see 
the caption of Table I 
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using HI of LDA band. 

The position of the /° peak appeared at about 3.2 eV below Ep. This is deep compared 
with the experimental result of about 2.4 eV. To obtain the f° level at about 2.7 eV, the 
occupation number of 4/ electron must be reduced to 0.8. This is very small compared with 
the value obtained by LDA calculation. We think that the total occupation of the 4/ electron 
should not so largely change from the value of LDA, because large electrostatic energy had 
been already included in the LDA band. The HI calculated directly from the band calculation 
seems to be somewhat large. 
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Fig. 10. Effective HI in DMFT for a-Ce. For the definitions of lines, see the caption of Fig. 4. 



5. CeSb 

CeSb has NaCl structure in which Ce ions form the fee lattice. The valence and the 
conduction bands mainly consist of the Sb 5p states and the Ce 5d states, respectively. They 
overlap slightly and CeSb has semi-metallic band structures. The top of the p valence band 
has the (j = 3/2)r§ character, and thus has the stronger hybridization matrix with the Tg 
CF state. 

In the p — f mixing model, a part of the valence band which has the same symmetry 
with the occupied 4/ band is pushed up above Ep due to the hybridization in the ordered 
phase. 31-34 Such re-constructed electronic band structures in the ordered states are expected 
by Hatree-Fock-like static mean field theory, 38 and confirmed by experimental studies. 35-37 
We need the DMFT band calculation in the paramagnetic phase. 

In Fig. 11, we show the HI for CeSb of LDA band, and in Fig. 12 the effective HI of DMFT 
calculation. The peak of HI in the energy region 0.25 ~ 0.3 Ry in Figs. 11 and 13 mainly 
corresponds to the j = 3/2 component of Sb p, and the peak at about 0.15 Ry corresponds 
to the j = 1/2 component. 26 ' 29 

In Fig. 11, we have peak structures near Ep which are caused by the hybridization of band 
states with the 4/ bands. The 4/ bands are pinned near Ep in the paramagnetic state. But 
the peaks do not necessarily appear just at Ep in the LDA band. We have sharp peaks just 
at Ep in the DMFT calculation as shown in Fig. 12. This may be caused by the hybridization 
with the correlated 4/ band. 

The occupation number ratio of the (j = 5/2)r7 to the (j = 5/2)r§, 0.39 is smaller than 
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Fig. 11. The HI calculated by the LDA band for CeSb. For the definitions of lines, see the caption 
of Fig. 4. 

ratio 1/2 of their degeneracy factor, thus the effective energy of the latter is lower in this case. 

In the single particle spectra shown by Fig. 13, a small hump is observed at about 0.7 eV 
below Ep. A large peak appears at about 3 eV below Ep. This double-peak-like structure is 
qualitatively similar to the result of experiments in PES. The peak structure at Ep in HI of 
DMFT causes the sharp peak just at Ep in the spectra. 

We show the magnetic excitation spectra by the bold dot-dashed line in Fig. 8. The 
excitation spectra for the fictitious model with matrix elements within the (j = 5/2)Fg man- 
ifold are shown by the thin dot-dashed line. The later decreases to the half of the low en- 
ergy value at the excitation energy —3.4 (63 K), while the former has contribution from the 
(j = 5/2)rg — > (j = 5/2)T7 excitation in the energy region about —3.3 (79 K). 

The mass enhancement factors are estimated to be very large. However, these large values 
will not have actual meaning as already noted because we have calculated them at the high 
temperature region, T = 300 K. 

At this point we should note that the present calculation for CeSb has started from the 
LDA band. Usually it leads to rather large DOS near Ep due to the large band overlapping 
of Sb p and Ce d bands. 38 ' 50 ' 51 Very high Tk will be partly ascribed to the use of this large 
DOS. In transport experiments, the Kondo temperature is expected to be about 10 K. 33 The 
correct estimation of the Tk will depend on this band overlapping. 

We tentatively shift the Sb 5s and bp bands —0.5 eV by putting wsb,p(— ) = <^Sb,s(— ) = 
-0.037 Ry . 51 The atomic energy level of 4/(j = 5/2)r 7 is also shifted -1000 K. This extra 
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Fig. 12. Effective hybridization in DMFT for CeSb. For the definitions of lines, see the caption of 
Fig. 4. 




Fig. 13. 4/ spectra of the auxiliary impurity model for CeSb. For the definitions of lines, see the 
caption of Fig. 5. 



shift may correspond to the CF energy splitting due to the electrostatic potential. The DMFT 
calculation for this modified model is carried out. The occupation numbers of (j = 5/2)Tj 
and (j = 5/2)r8 are obtained as 0.39 and 0.58, respectively. The effective CF energy becomes 
lower for (j = 5/2)Tj in this case. We note, however, that we obtain another type solution in 
which the occupied electrons are almost (j = 5/2)rs when the extra CF shift of (j = 5/2)Tj 
is reduced to —500 K. Both types of states are competing. 
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Fig. 14. 4/ spectra calculated by reducing the band overlapping between Sb p and Ce d for CeSb. The 
Sb p and s bands are shifted to low energy side 0.5 eV by putting Ldsb,s(— ) = ^Sb,s{— ) = —0.037 
Ry. The atomic energy level of 4f(j = 5/2)Tr is shifted —1000 K. For the definitions of lines, see 
the caption of Fig. 5. 



Table III. Various quantities obtained in DMFT calculation for CeSb. E-p = 0.3743 Ry, and energy 
levels are given in Ry. e\f — —0.16022 Ry. £4/ = 0.007043 Ry. For the definition of quantities, see 
the caption of Table I 





r 7 


r 8 


3 = 7/2 




0.3788 


0.4994 


0.3178 


„( im p'-) 


0.2753 


0.6917 


0.0458 


^( band ) 


0.2652 


0.6892 


0.0421 


4 imp ° 


0.2658 


0.6900 


0.00407 




-0.15815 


-0.15827 


-0.113914 




23.1 


29.9 


2.9 


£r 


0.17443 


0.13428 


0.28134 


f r 


0.1506 


0.1892 


0.0006 



The single particle spectra are shown in Fig. 14. They show the shallow energy peak at 
about 0.8 eV and the deep energy peak at about 3 eV below Ep. This result well corresponds 
to the experimental result. The magnetic excitation is shown in Fig. 8 by the bold two- 
dots-dashed line. The excitation of the fictitious model with the matrix elements within the 
(j = 5/2)T7 manifold is shown by the thin two-dots-dashed line. The later decreases to the half 
of the low energy value at about excitation energy —3.9 (20K), while the former has a peak 
at the excitation energy about —3.2 (100K) corresponding to the (j = 5/2)^ — ► (j = 5/2)Tg 
transition. We note that the peak of HI at E-p is considerably small and sharp compared with 
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the result in Fig. 12. It becomes very similar to the peak in HI used by Takeshige to reproduce 
T K of CeSb. 42 

Detailed systematic studies on PES and Fermi surface structures are necessary in quan- 
titative calculation. However, we note again, the hybridization with the correlated 4/ band 
will produce sharp peaks at Ep in HI. The similar result has been also found in ref. 43. The 
such effect has not been considered in previous studies. 26 It is also necessary to analyze the 
transport phenomena by including the contribution of the correlated 4/ band. 

6. Summary and Discussion 

We have developed a band calculation method based on the DMFT + LMTO frame work, 
and applied it to Ce metal and CeSb. The starting band is calculated by the LDA method, 
and the auxiliary impurity problem is solved by the NCA/ 2 v' method which is explained in 
Appendix. The method can include the f 2 configuration and gives the Kondo temperature 
(Tk) not the so different from that obtained by NRG, thus we can expect quantitative cal- 
culation. The spin-orbit (SO) splitting and the crystalline field (CF) splitting are included in 
the self-energy term. Tk is estimated from the calculation of the magnetic excitation spectra. 

In Ce metal, the total hybridization intensity (HI) near Ep is reduced in DMFT. Reflecting 
this the SO side peak in PES becomes relatively conspicuous compared to the calculation 
directly using HI in the LDA band. The difference of HI between (j = 5/2)Tj and (j = 5/2)Fg 
states leads to the energy level difference between them in the auxiliary impurity model. But 
in magnetic excitation, this difference does not directly appear because the HI itself causes 
the CF splitting. The Kondo temperature for 7-Ce is estimated to be about 70 K, and the 
CF splitting is smaller compared with this in the present calculation. This point should be 
checked carefully, however, because we have not considered the CF splitting energy due to 
the electrostatic potential. The Tk in a-Ce is estimated to be about 1700K. In experiments 
the magnetic excitation in a Ce is expected to be about 10 3 K. 49 The SO splitting in DMFT 
calculation for Ce metal is larger than Tk, thus the Kondo temperature is greatly reduced 
from that estimated by calculation neglecting the SO effect. The studies of a-j transition 
including the effect of SO and CF splittings are retained to the future. 

For CeSb, we have obtained similar results to results of the previous calculation based on 
the impurity model in the gross sense. The double-peak-like structure in PES experiments is 
reproduced. We have a peak very near Ep in HI. It leads to the higher Tk- The peak is caused 
by the hybridization with the correlated 4/ band in DMFT theory. 

When the band overlapping between the Sb p and Ce d is reduced, the double peak 
structure of PES and Tk about 20 K were obtained. In transport experiments, the Kondo 
temperature is expected to be about 10 K. 32 Systematic studies of band structure such as 
analysis of the dHvA effect, 35 ' 36 the angle-resolved-PES 40 and transport properties 32 will be 
carried out in the near future. 
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In both materials, the HI which is directly calculated by the LDA band calculation seems 
to be larger than the magnitude expected from the comparison with PES experiment. But in 
the DMFT calculation the effect of HI is relatively weakened. The calculated HI is somewhat 
lager, but not so drastically large than the magnitude expected from experimental data. 

The auxiliary impurity problem is solved by the NCA/ 2 v' method. It gives not so different 
excitation spectra from those calculated by NRG. However the NCA/ 2 v gives the improved 
result as shown in Fig. 2. The implementation of this method will also be desirable though it 
will need huge computation time. 
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Appendix: NCA/V method 

Let us define the Hubbard X operator between the atomic states a and 6, 

X(a,b) = \a >< b\. (A-l) 

We denote by |0 > the f° state, and by |n > the states in which the n orbital is occupied. 
When the two electrons occupy the r] and rj orbitals, the state is denoted as \r\r( >. The 
atomic state of the 4/ electron is expressed as 

H t = Y,E(a)X(a,a), (A-2) 

a 

where E(a) is the energy of the atomic state. The creation operator of the / electron is also 
expressed as 

/, + =<r,|/+|0>%0)| £ <Vv'\f^W>X({ VV '), V '). (A3) 
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Fig. A-l. Typical diagrams in the resolvent method, top: the second order diagram for the resolvent 
self-energies. The wavy line denotes the f° resolvent, the dashed line and the double dashed line 
denote the f 1 and the f 2 resolvents, respectively, middle: the fourth order diagrams for the self- 
energy of f° resolvent, bottom: higher order diagrams for the f° resolvent, which lead to the 
vertex correction. 

The resolvent for the atomic state a, R(a; z) is expressed by using the self energy of the 
resolvent T,(a; z) as 12,14 ' 17 

R(a- z) = — 1 -. (A-4) 

z — E(a) — L(a; z) 

The skeleton diagrams of the NCA approximation which include the f 2 configuration are 
given in the top panel of Fig. A-l. 14 The original NCA for f° — f 1 fluctuation model is given 
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by the diagram a and b. In the diagram a for the self-energy of the f° state, the fermion 
loop which consists of the solid and the dashed lines appears. It leads to the summation on 
the degeneracy factor, and has the large factor Nf. Based on this fact, the self-consistency 
solutions of a and b types are derived. A part of the fourth order terms for the / self-energy 
is given in the middle panel of Fig. A-l. The diagrams a and b can be included by the skeleton 
diagram a of the top panel. However, the diagram c cannot be included in it. This diagram 
has also two loops, so it will have Nj same as the a and b have. In the higher order diagrams, 
we must include the diagrams shown in the bottom panel of Fig. A-2. Diagrams of this type 
for the f° self-energy can be summed up when we solve the integral equation for the vertex 
part shown graphically by top or middle panels of Fig. A-2, 17,19 

A(r)- z, x) = 1 + J ' dx'M ri/ (x')R(i]i]';z + x + x')R(r]';z + x')f(x')A(r]';z,x'), (A-5) 

where we have used the type a equation in the Fig. A-2. The quantity M v (x) denotes the 
hybridization intensity v 2 p c (x). 
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Fig. A-2. Graphical representations of the integral equation for the vertex function (top and middle 
panel), and the equation for the resolvent self-energy of f° (bottom). 

The graphical expression of the /° self-energy is shown in the bottom panel in Fig. A-2, 
and the equation is given as, 17 

S(0; z)=J2 dxM v( x ) R (v, z + x)f(x)A(ri; z, x). (A-6) 
v J 

When one put A(r/; z, x) = 1 by neglecting the contribution from the f 2 state, eq. (A-6 ) 
corresponds to the equation in the usual NCA of (/° — f 1 ) fluctuation. 
The self-energies of the \i] > state and |r/r/ > are given 17 

+ J dxA(r];z - x,x)R(0;z - x)A(r];z - x,x)f(-x), (A-7) 
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xR(0; z-x- x?)A(rf; z-x-x', x')R(r}; z - x)f{-x)f{-x). (A-8) 



Fig. A-3. Graphical representation for the resolvent self-energies of J 1 (left panel) and f 2 (right panel) 
states. 

The graphical expression is shown in Fig. A-3. In ref. 17, it has been pointed out that 
the exchange coupling due to the virtual excitation to the f 2 can be led even when the x 
dependence of A(r/; z, x) is neglected. The equation for the vertex is simplified as 

A( V ; z) = 1 + R(2)* I *dx , M v (x')i?(r/';z + x')/(a;')A(r ? ';z). (A-9) 

where R(2)* is some representative value of the resolvent of the f 2 state for z ~ £f. We usually 
put as R(2)* = l/(e f - {2e f + 17)). The equations (A-6) ~ (A-9) and (A-4) are solved. As 
discussed in ref. 17, this scheme gives the exchange constant consistent with that of the S-W 
transformation in the leading order of the 1/Nf expansion. We note that the original NCA is 
correct to the next order of the 1/Nf expansion if the valence fluctuation is restricted within 
f l -f- 

The single particle excitation spectrum is given by calculating Green's function defined 

by 

G v = - < T T f v (T)f+ > . (A-10) 

By substituting eq. (A. 3) into this expression, it is decomposed to Green's function defined as 

G(ab, b'a'; r) = - < T T X(ab; T)X{b'a!) > . (A-ll) 

The diagonal part of the spectral function for the r\ particle excitation is given as the convo- 
lution integral 

9G(0r/, r/0 : x + 10) = 1 + ^ - J dy£(0; y)(<3R(v, V + x + i0)), (A-12) 

^G{rj[{rrrf), (w'W ■ x + iO) = 1 + ^ - J dy£(r/; y)(SSR(rrf;y + x + 10)), (A-13) 
where the quantities £(a;x) is called as the defect part defined as 

£(a; x) = -^e- px (--^R(a; x + iO)). (A-14) 

Zf IT 
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Fig. A-4. Graphical representation of the single particle excitation. The top panel is the NCA/ 2 v' 
part, and the bottom is off-diagonal part which are considered in NCA/ 2 v method. 

The graphical expression is given by the top panel in the Fig. A-4. We note that the 
spectrum intensity has the off-diagonal type term shown by the bottom in the Fig. A-4, which 



It is known that the integration by x of (A. 15) becomes zero. 17 This term transfers spectral 
intensity in the high energy region to the low energy region. The expressions (A. 12), (A. 13) 
and (A. 15) have convenient form to calculate the particle excitation part, i. e., x > 0. The 
expression can be transformed into a suitable form to calculate the hole excitation case (x < 
) by changing the integration variable from y to y + x. 14 



is given as 



^G(rf(rjrf),rfi; x + iO) = %G(0r}, {rpf)rj[; x + iO) 



vr^/ J J 
xQ(R(r]'; y + iO)A(r/; y - x + i0)i?(0; y - x' + iO)) 




x%(R(ri; y + x-x' + iO)R(r]r]'; y + x + iO)). 



(A-15) 
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